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Water  temperature  is  given  in  degrees  Celsius 
degrees  Fahrenheit  (°F)  by  the  following  equation: 


(°C),  which  can  be  converted  to 


°F  =  1.8(°C)+32. 

Selenium  concentrations  are  given  in  micrograms  per  liter  (yg/L) .  One  thousand 
micrograms  per  liter  is  equivalent  to  1  milligram  per  liter  (mg/L) .  Micrograms  per 
liter  is  equivalent  of  "parts  per  billion."  Isotope  composition  is  expressed  in 
parts  per  thousand  or  "per  mil." 
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Sea  level:  In  this  report,  "sea  level"  refers  to  the  National  Geodetic  Vertical 
Datum  of  1929  (NGVD  of  1929) — A  geodetic  datum  derived  from  a  general  adjustment  of 
the  first-order  level  nets  of  both  the  United  States  and  Canada,  formerly  called  Sea 
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ABSTRACT 

Salinity  and  selenivun  concentrations  in 
shallow  ground  water  of  the  western  San 
Joaquin  Valley,  California,  are  related 
to  the  geomorphology  and  hydrology  of  the 
alluvial  fans.  The  highest  salinity  and 
selenivun  concentrations  in  shallow  ground 
water  occur  in  alluvium  deposited  by 
ephemeral  streams  and  at  the  margins  of 
the  major  alluvial  fans,  where  there 
were  naturally  saline,  fine-grained 
soils.  Low-to-moderate  salinity  and 
selenium  concentrations  in  shallow  ground 
water   are   associated   with   upper   and 


middle  areas  of  the  major  alluvial  fans 
deposited  by  intermittent  streams.  Areas 
with  the  most  naturally  saline  soils  have 
been  irrigated  in  the  last  40  years. 
These  are  now  the  areas  of  highest  salin- 
ity and  selenium  concentrations  in  the 
shallow  ground  water.  Present-day  (1986) 
ground-water  salinity  is  significantly 
correlated  with  natural  soil  salinity  in 
these  areas.  Isotopic  data  indicate  that 
the  highest  salinity  and  selenium  con- 
centrations in  ground  water  present  at 
low  altitudes  resulted  from  evaporative 
concentration  when  the  water  table  was 
shallow. 
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INTRODUCTION 

Saline  soil  and  ground  water  have 
caused  water-management  problems  for 
agriculture  in  the  western  San  Joaquin 
Valley,  California,  for  decades.  Recent 
findings  of  high-selenium  concentrations 
in  agricultural  drain  water  and  waterfowl 
deformities  where  drain  water  was  dis- 
charged (Presser  and  Barnes,  1985; 
Deverel  et  al.  ,  1984,  Ohlendorf  et  al., 
1986)  have  increased  environmental  con- 
cern about  agricultural  water  management 
in  the  area.  The  high-selenium  drain 
water  comes  from  subsurface  drainage 
systems  that  remove  shallow  ground  water 
to  keep  the  water  table  low  enough  for 
crop  production.  Attention  is  being 
directed  toward  management  strategies  to 
decrease  drain-water  flows  and  selenium 
concentrations  in  ground  water.  Develop- 
ment of  effective  water-management  strat- 
egies for  the  western  valley  requires  an 
understanding  of  processes  affecting 
salinity  and  selenium  in  shallow  ground 
water. 

Previous  studies  dealt  mainly  with 
drained  areas  where  the  water  table  was 
generally  less  than  3  m  below  land  sur- 
face. Deverel  and  Millard  (1988)  and 
Deverel  et  al.  (1984)  evaluated  the  dis- 
tribution of  and  interrelations  among 
chemical  constituents  in  shallow  ground 
water  in  areas  of  the  western  San  Joaquin 
Valley  where  drainage  systems  have  been 
installed.  The  most  saline  and  selenium- 
enriched  ground  water  occurred  in  Coast 
Range  alluvial-fan  deposits,  which  in- 
cludes the  Little  Panoche  Creek,  Panoche 
Creek,  Cantua  Creek,  and  Los  Gatos  Creek 
alluvial  fans,  and  numerous  small  fans  in 
between.  Selenium  and  other  mobile  oxy- 
anions  (molybdenum,  vanadium,  and  boron) 
are  correlated  with  ground-water  salinity 
in  Coast  Range  alluvium  (Deverel  and 
Millard,  1988) .  Selenium  concentrations 
were  lower  and  were  not  correlated  with 
salinity  in  ground  water  underlying  the 
deposits  derived  mainly  from  Sierra 
Nevada  rocks. 


The  objective  of  this  study  was  to 
identify  the  geochemical,  hydrologic,  and 
geomorphic  processes  that  affect  the 
regional  distribution  of  salinity  and 
selenium  in  shallow  ground  water  of  the 
central  part  of  the  western  valley.  In 
this  report,  regional  distributions  of 
various  constituents  in  ground  water  are 
evaluated  in  relation  to  the  geomor- 
phology  and  hydrology  of  the  alluvial 
fans.  This  evaluation  serves  to  identify 
the  processes  that  affect  the  distribu- 
tion of  salinity  and  selenium.  Data 
presented  by  Deverel  et  al.  (1984)  and 
additional  data  collected  primarily  up- 
slope  of  areas  that  now  have  drainage 
problems  were  combined  and  analyzed.  The 
location  of  the  study  area  and  areal 
distribution  of  wells  sampled  for  chem- 
ical analyses  are  shown  in  figure  1. 
Additional  data  collected  from  28  wells 
installed  in  a  drained  agricultural  field 
(site  M2,  shown  in  fig.  1)  also  were 
analyzed. 

This  study  is  part  of  a  comprehensive 
investigation  of  the  hydrology  and  geo- 
chemistry of  the  western  San  Joaquin 
Valley  by  the  U.S.  Geological  Survey. 
The  studies  are  being  done  in  cooperation 
with  the  San  Joaquin  Valley  Drainage 
Program  and  as  part  of  the  Regional 
Aquifer  System  Analysis  Program  of  the 
U.S.  Geological  Survey. 


GEOHYDROLOGIC   SETTING 

California's  San  Joaquin  Valley  is  a 
trough  between  the  Sierra  Nevada  on  the 
east  and  the  folded  and  faulted  Coast 
Ranges  on  the  west.  The  Sierra  Nevada,  a 
fault  block  that  dips  southwestward,  is 
composed  of  igneous  and  metamorphic  rock, 
mostly  of  pre-Tertiary  age.  The  Diablo 
Range  of  the  California  Coast  Range, 
which  borders  the  study  area  to  the  west, 
consists  of  an  exposed  Cretaceous  and 
Upper  Jurassic  core  assemblage  of  marine 
origin  overlain  and  juxtaposed  with 
Cretaceous  marine  and  Tertiary  marine  and 
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FIGURE  1 .    Location  of  study  area  and  areal  distribution  of  wells  sampled  for  chemical  analyses. 


continental  deposits.  The  alluvial-fan 
deposits  of  the  western  valley  are 
derived  exclusively  from  Diablo  Range 
rocks  (Miller  et  al.,  1971).  The  basin 
deposits  are  derived  from  the  Sierra 
Nevada  and  the  Diablo  Range,  but  are 
dominated   by   Sierra   Nevada   deposits. 

The  alluvial  fans  on  the  westernmost 
edge  of  the  valley  were  deposited  by 
intermittent  and  ephemeral  streams  (Bull, 
1964a,  b)  .  Intermittent  streams  origi- 
nate in  the  mountains  of  the  Diablo  Range 
and  flow  for  several  months  each  year. 
Ephemeral  streams  originate  in  the  foot- 
hills of  the  Diablo  Range  and  flow  only 
during  storms.    The  fans  deposited  by 


ephemeral  streams  are  smaller  than  and 
encircled  by  the  four  major  alluvial  fans 
deposited  by  intermittent  streams — Little 
Panoche,  Panoche,  Cantua,  and  Los  Gatos 
Creeks  (fig.  1)  .  Shallow  deposits  of 
the  major  alluvial  fans  are  typically 
coarse  grained  in  the  upper-  and  middle- 
fan  areas  and  fine  grained  in  the  lower- 
fan  areas,  including  the  margins  of  the 
fans  (Harradine,  1950) .  The  fine-grained 
soils  at  the  margins  of  the  fans  gen- 
erally are  poorly  drained  and  often 
saline. 

Under  natural  conditions,  ground-water 
recharge  was  primarily  from  infiltration 
of  water  from  intermittent  streams  in  the 
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upper  parts  of  the  alluvial  fans  (Belitz, 
1988) .  Ground  water  flowed  from  south- 
west to  northeast,  reflecting  the  general 
topographic  trend  of  the  area.  Discharge 
from  the  system  was  by  evapotranspiration 
and  streamflow  along  the  valley  trough 
(Mendenhall  et  al.,  1916,  Davis  and 
Poland,  1957) .  Flowing  wells  mapped  by 
Mendenhall  et  al.  (1916)  indicate  that 
artesian  conditions  prevailed  along  a 
broad   stretch   of   the   valley   trough. 

The  predevelopment  ground-water-flow 
system  was  characterized  by  horizontal 
flow  from  the  recharge  area  to  the  valley 
trough.  Pumpage  of  ground  water  for 
irrigation,  which  began  in  the  late 
1800 's  and  reached  a  maximum  in  the 
1950' s,  substantially  altered  the  ground- 
water-flow  system.  Percolation  of  irri- 
gation water  past  crop  roots  throughout 
the  western  valley  became  the  primary 
mechanism  for  ground-water  recharge. 
Evapotranspiration  by  crops  and  discharge 
from  wells  became  the  primary  discharge 
mechanisms.  The  combination  of  ground- 
water withdrawals  and  increased  recharge 
resulted  in  vertical  gradients  in  the 
regional  ground-water  system,  causing  a 
substantial  component  of  downward  flow. 

Ground-water  withdrawals  until  1952  in 
the  area  north  of  Mendota  and  until  1967 
in  the  area  south  of  Mendota  caused  the 
water  table  to  decline  along  the  fan 
margins  and  in  the  valley  trough.  When 
surface  water  imported  from  northern 
California  replaced  ground  water  for 
irrigation,  the  water  table  rose  through- 
out the  area.  The  percentage  of  the  area 
characterized  by  a  water  table  within 
6  m  of  land  surface  has  increased 
substantially  since  1952  (Belitz,  1988) . 


METHODS 

Differences  in  chemical  compositions  of 
ground  water  were  evaluated  using  Piper 
diagrams  (Hem,  1985)  and  the  chemical 
equilibrium  model  WATEQF  (Plummer  et  al. , 
1976) .  Mineral  phases  affecting  the 
ground-water  chemistry  were  assessed 
using  mineral  saturation  indices  deter- 
mined by  the  WATEQF  chemical  equilibrium 


model,  where  the  saturation  index  is 
equal  to  the  quotient  of  the  ion-activity 
product  and  mineral-solubility  constant. 
Negative  saturation  index  values  indicate 
mineral  undersaturation  and  positive 
values  indicate  supersaturation. 

The  hydrologic  history  of  selected 
water  samples  was  assessed  from  the 
deterinination  of  tritium  concentrations 
and  deuterium  and  oxygen-18  enrichment. 
Tritium  concentrations  are  expressed  in 
tritium  units  (TU)  ,  where  1  TU  is  equal 


to  one 
atoms . 


tritium  atom  in  10 


li 


hydrogen 


The  distribution  of  various  chemical 
constituents  in  the  ground  water  was 
compared  using  geostatistical  methods. 
Kriging  was  used  to  generate  interpolated 
values  of  three  constituents  (ground- 
water salinity,  soil  salinity,  and  sele- 
nium in  ground  water)  based  on  the  spa- 
tial variance  structure  and  the  measured 
values.  The  estimated  variance  structure 
is  defined  by  the  semivariogram  (Journel 
and  Huijbregts,  1978) : 


Y(h)  =  iE[(z(x  )-z(x.))2], 
i     3 


(1) 


where 


h  = 


x.-x.  is  the  distance  (lag)  between 
pairs  of  data, 

z  is  the  logjQ  value  of  the 
constituent  at  site  x.  or 
x  . ,  and 

: 

E  designates  the  expected  value 
of  the  terms  inside  the 
brackets. 


The  frequency  distributions  of  the 
constituents  evaluated  (ground-water 
salinity,  expressed  as  dissolved  solids; 
soil  salinity;  and  selenium  in  ground 
water)  were  log  normal.  This  was  deter- 
mined by  comparing  the  log  normal  scores 
to  the  logarithms  (base  10)  of  the  values 
of  the  constituents.  The  log  normal 
scores  were  significantly  correlated  with 
the  logarithms  of  the  values.  The  cor- 
relation coefficients,  r,  for  dissolved 
solids  in  ground  water,  soil  salinity, 
and  selenium  in  ground  water  are  0.97, 
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0.98,  and  0.98,  respectively.  All  r 
values  are  significant  at  a=0.01.  The 
logarithms  of  the  data  are  used  in  the 
geostatistical  analysis  (Delhomme,  1978; 
Neuman,  1982) .  The  semivariogram  was 
verified  by  sequentially  estimating  each 
measured  value  with  the  remaining  data. 
Kriging  of  historical  soil  salinity  in 
the  western  valley  was  done  with  573 
point  estimates  of  subsurface  soil 
salinity  from  Harradine  (1950) . 


Field 


Samples  of  shallow  ground  water  were 
collected  at  118  observation  wells  dis- 
tributed throughout  the  study  area  (fig. 
1)  .  Seventy-two  of  the  wells  were  sam- 
pled during  an  18-day  period  in  spring 
1984.  Forty-six  wells  were  sampled  in 
the  spring  of  1986.  Samples  collected  at 
these  46  wells  were  analyzed  for  deuter- 
ium and  oxygen-18  (fig.  1)  .  The  wells 
are  6  to  10  m  deep  and  constructed  of 
polyvinyl  chloride  pipe  perforated  over 
the  entire  depth  of  the  well. 

An  additional  28  wells  with  30-cm 
screens  at  3,  6,  10,  and  12  m  below  land 
surface  were  installed  at  the  M2  site 
(fig.  1) .  The  wells  were  installed  along 
a  300-m  transect  in  January  1986  and  were 
sampled  in  spring  1986. 

Most  of  the  observation  wells  were 
sampled  using  tubing  and  positive- 
displacement  bladder  pumps  constructed  of 
Teflon^.  The  submersible,  2-cm  diameter 
bladder  pump  operates  by  allowing  water 
to  enter  through  a  one-way  valve  in  the 
bottom  of  the  pvimp.  Pressurized  nitrogen 
gas  is  cycled  in  and  out  of  a  Teflon 
bladder  inside  the  pump  casing,  dis- 
placing the  water  up  the  tubing.  Some  of 
the  observation  wells  were  sampled  using 
a  peristaltic  pump  to  remove  water  from 
the  well  through  Teflon  tubing  into 
sample  containers. 


Prior  to  sampling,  each  well  was  pumped 
until  the  specific  conductance  of  the 
water  did  not  vary  by  more  than  10  per- 
cent for  three  consecutive  well-casing 
volumes  of  water.  Electrical  conduc- 
tivity meters  were  standardized  using 
potassium  chloride  standards  within  2,000 
US/cm  of  the  well  water.  Specific  con- 
ductance and  pH  were  measured  several 
times  at  each  well  prior  to  collection  of 
samples.  Before  each  sample  was  taken, 
another  well-casing  volume  was  pumped 
with  the  bladder  pump  only. 

Prior  to  sampling  at  each  site,  all 
sampling  apparatus  and  containers  were 
rinsed  three  times  with  the  water  to  be 
sampled.  Well  water  was  used  to  rinse 
filters,  filter  holders,  and  sample  con- 
tainers. Prior  to  collection  of  samples 
for  laboratory  analysis  of  the  inorganic 
constituents,  unfiltered  50-mL  samples 
were  titrated  incrementally  with  dilute 
sulfuric  acid  for  bicarbonate  and  car- 
bonate concentrations.  The  pH  was  mea- 
sured with  portable  meters  that  were 
calibrated  at  each  site  with  standardized 
solutions.  Portable  pH  meters  and  elec- 
trodes were  calibrated  at  each  well  and 
for  each  bicarbonate  titration.  The 
calibrations  were  done  using  pH  buffers 
maintained  at  the  temperature  of  the 
sample  water  within  ±2  °C.  Temperature 
was  measured  with  mercury  thermometers 
that  were  checked  against  a  standardized 
laboratory  thermometer. 

Samples  for  analysis  of  dissolved  con- 
stituents were  pressure  filtered  through 
0.45-tim  membrane  filters.  Samples  for 
determination  of  major  ions  and  all  trace 
elements  were  collected  and  stored  in 
polyethylene  bottles  with  Teflon-lined 
caps.  All  samples  were  acidified  with 
nitric  acid  to  a  pH  of  less  than  2  except 
those  collected  for  determination  of 
chloride,  sulfate,  bicarbonate,  carbon- 
ate, pH,  and  specific  conductance,  which 


^The  use  of  brand  names  in  this  report  is  for  identification  purposes  only 
and  does  not  constitute  endorsement  by  the  U.S.  Geological  Survey. 


Methods 


were  untreated.  Samples  collected  for 
the  determination  of  selenate  and  sele- 
nite  were  acidified  to  a  pH  less  than  1 
using   concentrated   hydrochloric   acid. 

Determination  of  the  oxidation- 
reduction  status  of  the  ground-water 
samples  was  done  by  determining  the 
dissolved-oxygen  concentration,  the 
platinum  electrode  redox  potential,  and 
the  ratio  of  selenite  to  selenate  acti- 
vities. The  oxidation-reduction  poten- 
tial was  calculated  using  the  Nernst 
equation  and  Eo  values  (Bard  et  al., 
1985)  for  the  reactions 


0^  +  4H    +  4e     -••  2H2O, 


SeO.    ■<-  +  2H     +  2e     *  SeO,  +  +  H,0. 

H    ■>■  +3+2 


(2) 
(3) 


The  platinum  electrode  redox  potential 
and  dissolved  oxygen  were  determined  by 
placing  the  platinum  electrodes  in  a  flow 
through  chamber  as  described  by  Wood 
(1981) .  Samples  for  analysis  of  oxygen- 
and  hydrogen-stable  isotopes  and  tritium 
were  unfiltered  and  collected  in  glass 
bottles  with  polyethylene  lids  so  that 
there  was  no  air  space . 


digestion,  together  with  the  inorganic 
selenium  originally  present,  is  then 
reduced  to  the  selenite  form  using  a 
stannous  chloride  and  potassium  iodide 
mixture.  The  selenium  hydride  is  gen- 
erated by  reducing  the  selenite  form 
using  sodium  borohydride.  The  hydride 
gas  is  stripped  from  the  solution  by  a 
stream  of  nitrogen  gas  and  its  concentra- 
tion is  determined  by  atomic-absorption 
spectrometry . 

To  determine  the  concentrations  of 
selenate  and  selenite,  three  separate 
analyses  were  done  on  each  sample. 
First,  the  total  concentration  of  the 
dissolved  forms  was  determined  by  the 
above  method.  Second,  the  sample  was 
analyzed  without  performing  the  oxidative 
digestion.  Third,  the  sample  was 
analyzed  without  performing  the  stannous 
chloride  reduction.  The  second  analysis 
determines  the  concentration  of  the  in- 
organic selenium  and  the  third  analysis 
determines  the  concentration  of  the 
selenate  species.  The  selenite  concen- 
tration was  calculated  from  the  differ- 
ence between  the  results  of  the  second 
and  third  steps. 


Laboratory 

Calciimi,  magnesium,  sodium,  potassium, 
sulfate,  chloride,  aluminum,  iron,  and 
silica  were  determined  by  atomic- 
absorption  spectrometric  methods  des- 
cribed in  Skougstad  et  al.  (1979)  . 
Chloride  was  determined  by  the  ferric 
thiocyanate  colorimetric  method.  Sulfate 
was  determined  by  the  turbidimetric 
procedure. 

Selenium  was  determined  by  hydride 
generation  and  atomic  absorption  spec- 
trometry (Fishman  and  Bradford,  1982)  . 
The  method  used  for  selenium  in  this 
study  is  designed  to  determine  the  total 
concentration  of  all  forms  of  selenium 
present  in  the  water  sample.  A  sample  is 
first  subjected  to  an  oxidative  digestion 
to  release  any  selenium  from  the  organic 
fraction.   The  selenium  released  by  this 


The  oxygen-isotopic  compositions  of  the 
water  samples  were  determined  using  the 
carbon  dioxide  equilibration  method 
(Epstein  and  Mayeda,  1953).  Analysis  was 
done  with  a  60°  sector,  150-mm  radius 
double-collecting  isotope  ratio  mass 
spectrometer.  The  results  are  reported 
relative  to  Vienna  Standard  Mean  Ocean 
Water  (V-SMOW) : 


180/160        -  180/160 

sample V-SMOH   ,  „„„  ,,. 

6180  =  *^ X  1,000.  (4) 


(180/160) 


V-SMOW 


The  hydrogen-isotopic  compositions  were 
determined  by  analyzing  hydrogen  quanti- 
tatively extracted  from  water  (Kendall 
and  Coplen,  1985) .  Hydrogen  results  are 
reported  relative  to  V-SMOW  in  the  per 
mil  notation.  The  standard  deviation  of 
the  laboratory  measurements  of  oxygen- 
and  hydrogen-isotopic   compositions   are 
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0.10  and  1  per  mil,  respectively.  The 
hydrogen-  and  oxygen-isotopic  composi- 
tions were  determined  in  the  Isotope 
Fractionation  Project  Laboratory  of  the 
U.S.  Geological  Survey,  Reston,  Virginia. 

Tritixom  was  determined  by  scintillation 
counting  after  electrolytic  enrichment 
(International  Atomic  Energy  Agency, 
1976) .  All  tritium  analyses  were  done  by 
the  Environmental  Isotope  Laboratory  at 
the  University  of  Waterloo,  Ontario, 
Canada. 


RESULTS  OF  GEOSTATISTICAL 
ANALYSIS 

The  kriged  values  for  ground-water 
salinity,  soil  salinity,  and  selenium  in 
ground  water  were  used  to  evaluate  the 
relation  of  salinity  and  selenium  in 
ground  water  to  hydrologic  and  geochem- 
ical  processes.  The  theoretical,  iso- 
tropic, semivariograms  for  these  three 
constituents  were  estimated  with  spheri- 
cal models  (fig.  2)  .  The  values  of  the 
variables  in  the  equations  describing  the 
variograms  are  shown  in  table  1.  The 
spherical  equation  has  the  form 


■If(h)  =  iD[3/2(h/a)+i(h/a)3]+c, 


(5) 


where  co+c,  a,  and  c  are  the  sill,  range 
of  influence,  and  nugget  variance,  res- 
pectively (Journel  and  Huijbregts,  1978) . 
For  all  three  constituents,  the  semivar- 
iogram  was  calculated  in  all  directions 
(isotropic,  east-west,  north-south, 
northeast-southeast,  and  southwest- 
northeast)  .  No  drift  was  apparent  in  any 
direction  so  the  isotropic  semivariogram 
was  used  to  estimate  the  spherical  model 
to  describe  the  spatial  variance 
structure. 

The  semivariogram  model  was  verified  by 
estimating  the  data  for  each  of  the  three 
constituents  listed  above.   Each  of  the 


observations  was  removed  sequentially  and 
was  kriged  using  the  remaining  observa- 
tions. The  semivariogram  seems  reason- 
able because  the  mean  estimation  error 
(M  )  for  the  three  constituents  was  close 
to  zero  (table  1)  and  estimation  errors 
appear  normally  distributed  around  zero 
for  the  range  of  estimated  values.  Also, 
the  standard  deviation  of  the  reduced 
error  (o  )  is  close  to  one  (table  1)  . 
The  origxnal  observations  estimated  by 
this  procedure  were  significantly  corre- 
lated (a=0.01)  with  the  estimated  values 
for  ground-water  salinity,  soil  salinity, 
and  selenium  in  ground  water. 


The  estimation  error  (M  )  is  defined  as 

e 


M  =  E(z  -z), 
e     1 


(6) 


and  the  standard  deviation  of  the  reduced 

error  (a   )  is  defined  as 
RE 


E[(z^ 


-z,)VsJ^J 


where  E  denotes  the  expected  value, 
symbols  z.   ,z. 


RE      i   i    i]',         (7) 

The 
and  sf  represent  the 
estimated  value  of  the  logarithm  of  the 
constituent  (ground-water  salinity,  soil 
salinity,  and  selenium  in  ground  water) , 
the  actual  value  of  the  logarithm  of  the 
constituent,  and  the  variance  of  the 
estimation   at   point   i,   respectively. 

The  normality  of  the  estimation  error 
for  the  kriged  values  was  tested  by 
calculating  the  normal  scores  and  com- 
paring them  to  the  estimation  error 
values.  For  all  three  constituents,  the 
r  values  for  the  regression  of  the  esti- 
mation error  values  and  the  normal  scores 
were  significant  at  a=0.01.  The  r  values 
were  0.99  for  ground-water  salinity,  0.99 
for  soil  salinity,  and  0.96  for  selenixim 
in  ground  water.  The  normalized  standard 
deviation  of  the  kriged  estimate  (David, 
1977)  varied  normally  around  zero  for  the 
range  of  the  estimated  values  for  all 
three  constituents. 
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FIGURE  2.  Semivariograms  for  the  log^Q  of  (a)  dissolved  solids  in  ground 
water;  (b)  percent  soil  salts  from  Harradine(  1950);  and  (c)  selenium  in  ground 
water.  The  solid  lines  represent  the  spherical  estimates  of  the  semivariograms. 
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TABLE  l.—Semivariogram  model  variables,   mean  estimation  error,   and  standard 
deviation  of  the  reduced  error  for  logjo  of  ground-water  salinity,   soil 
salinity,   and  selenium  in  ground  water 

[Symbols   a,    c,    and  w  are   the  parameters   in  the 
equation  Y  (h)=(jo  [3/2  (h/a) +i  (h/a)  3]+c] 


Constituent 


Standard 

deviation 

of  the 

Mean 

Number  of 

reduced 

estimation 

data  points 

error 

error 

'\e^ 

(M  ) 
e 

111 

1.023 

0.0023 

Dissolved   solids   in 
ground  water    (mg/L) 

Surface   soil   salinity 
(percent) 

Selenium  in  ground 
water    (pg/L) 


573 


118 


,97 


.97 


.0026 


,0045 


15  0.07      0.21 


17.5         .06         .095 


.2  .46 


RELATION   OF   SALINITY  AND 

SELENIUM   TO   HYDROLOGIC   AND 

GEOCHEMICAL   PROCESSES 

Alluvial-Fan  Geomorphology  and 
Hydrology 

The  distribution  of  salinity  (expressed 
as  dissolved  solids)  in  shallow  ground 
water  reflects  the  geomorphology  and  hy- 
drology of  the  alluvial  fans,  which  were 
deposited  by  streams  originating  in  the 
Coast  Range  and  foothills.  Dissolved- 
solids  concentrations  generally  are  simi- 
lar for  the  three  largest  alluvial  fans 
associated  with  intermittent  streams.  In 
ground  water  in  the  Panoche  Creek  allu- 
vial fan,  dissolved  solids  generally  are 
less  than  3,000  mg/L  in  the  middle-fan 
areas  (fig.  3)  .  In  the  lower- fan  area, 
dissolved  solids  are  commonly  greater 
than  5,000  mg/L,  with  some  areas  at  the 


fan  margins  that  have  dissolved-solids 
concentrations  greater  than  10,000  mg/L. 
Dissolved  solids  also  generally  are  less 
than  5,000  mg/L  in  the  middle-fan  areas 
of  the  Cantua  and  Los  Gatos  Creek  allu- 
vial fans  and  generally  are  greater  than 
5,000  mg/L,  and  are  as  much  as  to  10,000 
mg/L  at  the  margins  of  the  fans.  Sparse 
data  prevent  evaluation  of  similar 
patterns  for  Little  Panoche  Creek  fan. 

Dissolved-solids  concentrations  in 
areas  of  ephemeral-stream  alluvial  fans 
generally  are  higher  than  at  similar 
positions  in  adjacent  intermittent-stream 
alluvial  fans.  Ground  water  in  the 
middle-fan  areas  of  the  ephemeral-stream 
deposits  north  of  Cantua  Creek  has 
dissolved-solids  concentrations  that 
generally  are  between  3,000  and  5,000 
mg/L  and  are  as  much  as  10,000  mg/L.  In 
the  lowest  altitudes  of  these  ephemeral- 
stream  deposits,  dissolved  solids  are 
commonly  greater  than  10,000  mg/L. 
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FIGURE  3.    Dissolved  solids  in  shallow  ground  water,  1986. 


Over  geologic  time,  smaller  quantities 
of  water  flowed  in  the  ephemeral  streams 
originating  in  the  foothills  than  in  the 
larger  intermittent  streams.  The  soils 
of  all  parts  of  the  alluvial  fans  assoc- 
iated with  these  ephemeral  streams  have 
been  subject  to  less  natural  leaching 
than  the  soils  of  the  middle  and  upper 
parts  of  the  alluvial  fans  of  the  inter- 
mittent streams.  Thus  the  soils  of  the 
ephemeral-stream  fans  were  historically 
more  saline  than  the  middle  and  upper 
parts  of  the  intermittent-stream  fans. 
Harradine  (1950)  estimated  percent  soil 
salts  in  the  1940's  at  573  points  that 
are  shown  in  figure  4.   Subsurface  soil 


salinity  based  on  Harradine  (1950)  and 
the  area  of  flowing  artesian  wells  de- 
fined by  Mendenhall  et  al.  (1916)  are 
shown  in  figure  5.  Soils  that  have  0.5 
percent  salt  or  greater  are  present  in 
the  middle-fan  areas  of  the  alluvial  fans 
that  were  deposited  by  ephemeral  streams. 
Saline  soils  that  have  greater  than  0.75 
percent  salt  extend  into  the  middle-fan 
areas  in  the  ephemeral-stream  fan  depos- 
its to  the  north  of  the  Cantua  Creek 
alluvial  fan.  In  contrast,  soils  of  the 
alluvial  fans  deposited  by  intermittent 
streams  are  less  saline  in  the  upper-  and 
middle-fan  areas,  with  soils  that  have 
0.5  percent  salt  or  less. 
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FIGURE  4.    Distribution  of  soil-sampling  sites. 


Low-lying  areas  of  high  salinity  and 
selenium  were  natural  ground-water- 
discharge  areas  (Mendenhall  et  al.,  1916; 
Davis  and  Poland,  1957) .  Figure  5  shows 
the  area  of  flowing  artesian  wells  des- 
cribed by  Mendenhall  et  al.  (1916)  in 
relation  to  soil  salinity  mapped  by 
Harradine  (1950)  .  The  proximity  and 
overlap  of  areas  of  highest  soil  salinity 
and  flowing  artesian  wells  indicates  that 
the  probable  reason  for  naturally  saline 
soils  was  evaporative  ground-water  dis- 
charge at  the  margins  of  the  alluvial 
fans. 

Salinity  of  ground  water  and  soils  also 
reflects  irrigation  patterns  since  the 


early  1900 's.  Figure  6  shows  the  areas 
that  were  irrigated  with  ground  water  as 
of  1940  (Bull  and  Miller,  1975) .  At  the 
time  of  Harradine 's  sampling  in  the 
1940 's,  the  areas  of  highest  salinity 
were  at  the  southern  and  northern  margins 
of  the  Panoche  Creek  fan,  and  at  the 
margins  of  the  Cantua  Creek  and  Los  Gatos 
Creek  fans  (fig.  5) .  These  saline  areas, 
which  include  the  middle-fan  areas  of  the 
ephemeral-stream  fans,  probably  were 
avoided  for  irrigated  agriculture  until 
after  1940  because  of  the  high  soil 
salinity.  Unirrigated  areas  at  higher 
altitudes  probably  were  not  irrigated 
because  of  greater  depths  to  ground 
water. 
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FIGURE  5.    Subsurface  soil  salinity,  1950. 


In  the  area  irrigated  after  1940,  the 
present-day  ground-water  salinity  is 
partly  the  result  of  leaching  of  soluble 
soil  salts  by  irrigation.  Figure  7  shows 
the  comparison  of  989  kriged  values  for 
log  10  of  ground-water  salinity  and  soil 
salinity.  Ground-water  and  soil  salinity 
are  significantly  correlated  (r=0.5, 
a=0.01) . 

The  texture  of  the  alluvial  deposits 
also  affects  the  relation  of  ground-water 
salinity  to  soil  salinity  in  this  area. 
About  50  estimated  points  denoted  in 
figure  7  are  in  a  single  area  south  of 
the  Cantua  Creek  alluvial  fan,  where 
coarse-textured  deposits  predominate 
below  3  m  (Julie  Laudon,  U.S.  Geological 
Survey,   written   commun. ,   1988) .     By 


removing  these  points,  which  generally 
represent  low  ground-water  salinity 
values  and  high  soil  salinity  values,  the 
resulting  correlation  coefficient,  r,  is 
0.66.  The  coarser  texture  of  the  depos- 
its of  this  area  relative  to  the  other 
deposits  in  the  area  irrigated  after  1940 
probably  resulted  in  more  rapid  leaching 
of  soluble  salts  by  irrigation  water  that 
percolated  through  the  root  zone. 
Present-day  ground-water  salinity  is 
uncorrelated  with  soil  salinity  as 
reported  by  Harradine  (1950)  (r=0.19)  in 
the  area  irrigated  before  1940.  The 
Wilcoxon  rank  sum  test  (Hollander  and 
Wolfe,  1973)  indicated  that  soil  salinity 
was  significantly  (a=0.01)  greater  in  the 
area  irrigated  after  1940. 


12   Salinity  and  Selenium,  Western  San  Joaquin  Valley,  California 


Kettleman 
City 


FIGURE  6.   Areas  irrigated  as  of  1940. 


Saline  soils  are  shown  in  an  area  irri- 
gated by  surface  water  prior  to  1940  at 
the  northern  margin  of  the  Panoche  Creek 
fan,  where  ground-water  salinity  is  now 
low.  Because  this  area  has  been  irri- 
gated with  surface  water  since  the  early 
1900 's,  an  irrigation-induced  rise  in  the 
water  table  probably  led  to  saline  con- 
ditions in  the  unsaturated  zone  because 
of  evapotranspiration.  The  soils  of  this 
area  were  mapped  as  nonsaline  in  1915 
(Nelson  et  al.,  1919).  Installation  of 
drainage  systems  in  this  area  since  1950 
has  allowed  displacement  of  saline  ground 
water  and  leaching  of  soil  salts.  This 
is  probably  the  reason  that  the  present- 
day  ground-water  salinity  of  the  area  is 
low. 


The  relative  distribution  of  selenium 
in  ground  water  is  similar  to  the  salin- 
ity distribution  (fig.  8)  ,  and  selenium 
in  ground  water  and  salinity  are  signif- 
icantly correlated  (Deverel  and  Millard, 
1988) .  Selenium  concentrations  generally 
are  less  than  20  yg/L  in  ground  water  in 
the  middle-fan  areas  of  the  alluvial  fans 
deposited  by  intermittent  streams.  In 
the  lower-fan  areas,  principally  at  the 
northern  and  southern  margins  of  the 
Panoche  Creek  fan,  concentrations  are  as 
much  as  several  hundred  micrograms  per 
liter.  Selenium  concentrations  in  ground 
water  in  the  middle-fan  areas  of  alluvial 
fans  deposited  by  ephemeral  streams  range 
from  50  to  500  ug/L.  High  selenium  in 
ground  water  generally  corresponds  to  the 
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FIGURE  7.   Comparison  of  estimated  log  jg  of  ground-water  salinity  and  soil  salinity 

for  the  area  irrigated  after  1940. 


saline  soil  areas.  Therefore,  leaching 
of  salts  from  saline  soils,  particularly 
in  areas  that  were  irrigated  after  1940, 
probably  contributed  substantial  amounts 
of  selenium  to  the  shallow  ground  water. 


Tidball  et  al.  (1986)  assessed  the 
distribution  of  selenium  in  soil  in  the 
western  San  Joaquin  Valley.  The  highest 
selenium  concentrations  (99th  percentile) 
were  in  the  upper-fan  areas  of  the 
ephemeral-stream  fans  between  the  Cantua 
Creek  and  Panoche  Creek  alluvial  fans. 
These  are  areas  that  have  been  irrigated 
only  since  the  1960 's  and  were  subjected 
to  little  natural  leaching.  High  soil 
selenium  concentrations  (95-99th  percen- 
tile) also  were  along  the  margin  of  the 
Panoche  Creek  fan  and  at  the  base  of  the 
ephemeral-stream  deposits  south  of  the 
Panoche  Creek  alluvial  fan.  These  high 
soil  concentrations  are  mostly  the  result 
of  evapotranspiration  from  a  shallow 
water  table. 


Studies  by  Amundson  et  al.  (1986)  and 
Fujii  et  al.  (1988)  provide  information 
about  the  nature  and  fate  of  the  total 
soil  selenium.  Amundson  et  al.  (1986) 
examined  changes  in  potassium  phosphate 
{K2HP0£j)  extractable  selenium  in  soils  on 
the  Panoche  Creek  alluvial  fan  from  1946 
and  1985.  Soil  samples  collected  from 
the  same  locations  in  1946  and  1985 
showed  substantial  decreases  in  K2HP0j^ 
extractable  selenium.  Soils  irrigated 
after  1940  were  higher  in  K2HP0ij  extract- 
able  selenium  in  1946  than  the  soil  irri- 
gated before  1940.  All  1985  K2HP0t^ 
extractable  selenium  values  were  about  10 
percent  of  the  1946  value.  Fujii  et  al. 
(1988)  also  found  that  selenium  in  irri- 
gated soils  is  presently  in  a  form 
resistant  to  leaching.  They  showed  that 
only  a  small  part  of  the  total  selenium 
in  soils  sampled  in  1985  in  fields  irri- 
gated for  more  than  30  years  is  in  water 
soluble  or  extractable  forms.  The  most 
mobile  forms  of  selenium  apparently  have 
been  mostly  removed  from  the  soils  in  the 
study  area  by  irrigation. 
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FIGURE  8.    Selenium  concentrations  in  ground  water,  1986. 


Evaporative  Concentration 

Although  the  leaching  of  natural  soil 
salts  by  irrigation  water  was  probably 
the  primary  process  that  initially  formed 
saline  shallow  ground  water,  subsequent 
evaporative  concentration  caused  further 
increases  in  salt  and  selenium  concentra- 
tions. Figure  9  shows  the  isotopic  com- 
position for  46  selected  samples  col- 
lected throughout  the  study  area 
(fig.  1)  .  Samples  collected  at  the 
lowest  altitudes  (50  to  58  m  above  sea 
level)  are  mostly  enriched  in  the  heavy 
isotopes  and  follow  an  evaporative  trend 
line.  Deverel  and  Fujii  (1988)  reported 
that  ground-water  isotopic  data  collected 
in  drained  agricultural  fields  fit  a  sim- 
ilar evaporative  trend  line.   At  higher 


altitudes,  where  selenium  concentrations 
generally  are  low  (fig.  8)  ,  the  ground 
water  is  more  depleted  and  has  been  sub- 
ject to  little  or  no  evaporation. 

Most  of  the  shallow  ground-water  sam- 
ples collected  at  the  higher  altitudes 
(71  to  100  m  above  sea  level)  represent  a 
local  meteoric  water  line  (Davis  and 
Coplen,  1988) .  The  isotopic  composition 
of  ground  water  collected  from  deep  parts 
of  the  regional  unconfined  aquifer  (at 
100  to  200  m  below  land  surface)  in  areas 
where  the  most  isotopically  enriched 
shallow  ground-water  samples  were  col- 
lected, is  similar  to  the  approximate 
intercept  of  the  evaporative  trend  line 
with  the  apparent  local  meteroic  water 
line  (fig.  9)  .  This  ground  water  was 
applied  as  irrigation  water  throughout 
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FIGURE  9.  Comparison  of  delta  deuterium  and  delta  oxygen-18 
for  samples  collected  from  various  altitudes. 


the  study  area  prior  to  1970  when  the 
water  table  was  too  deep  for  evaporative 
concentration  to  occur.  The  saline 
shallow  water  collected  at  the  lower 
altitudes  seems  to  be  regional  ground 
water  that  was  applied  as  irrigation 
water  and  then  evaporated  where  the  water 
table  was  shallow. 

The  two  samples  collected  at  low  alti- 
tudes had  tritium  concentrations  of  1.4 
and  1.7  TU  and  isotopic  compositions  of 
6d=-52.5  and  6lSo=-4.30  and  6d=-53.0, 
6^°O=-5.60  per  mil,  respectively.  The 
low  tritium  values  in  these  two  samples 
indicate  that  this  partially  evaporated 
water  is  composed  primarily  of  water  that 
infiltrated  prior  to  the  introduction  of 
atmospheric  tritium.  In  this  area,  where 
there  is  probably  little  ground-water 
recharge  from  rainfall,  we  assume  that 
most  of  the  tritiiim  in  the  ground  water 
was  introduced  since  about  1970,  when 
water  from  the  Federal  Central  Valley 
Project  was  first  introduced.  The 
tritium  values  in  13  samples  from  higher 


altitudes  ranged  from  7.5  to  31.5  TU 
(median,  15.25  TU)  ,  indicating  that  the 
ground  water  infiltrated  at  land  surface 
since  about  1970.  Irrigation  water  col- 
lected from  the  Federal  Central  Valley 
Project  in  spring  1986  had  13  TU. 

The  most  probable  explanation  for  the 
lower  tritium  values  at  the  lower  alti- 
tudes is  the  greater  depth  below  the 
water  table  at  which  the  samples  were 
collected.  All  wells  sampled  are  perfo- 
rated over  their  entire  length  such  that 
wells  at  the  lower  altitudes  penetrate 
farther  into  the  saturated  zone.  Wells 
sampled  at  the  lower  altitudes  probably 
intercept  older  evaporated  water.  Wells 
at  the  higher  altitudes  intercept  only 
the  top  1  to  3  m  of  the  ground  water. 

Clusters  of  wells  that  were  installed 
at  the  M2  site  (fig.  10)  at  various 
depths  and  screened  over  small  intervals 
(0.3  m)  provide  additional  information 
about  the  processes  of  leaching  and  evap- 
oration.   Figure  10  shows  a  comparison 
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FIGURE  10.   Comparison  of  delta  deuterium  and  delta  oxygen-18 
for  samples  collected  at  various  depths  at  the  M2  site. 


of  i5d  and  6^^0  for  samples  collected  from 
wells  along  a  300-m  transect  in  the  agri- 
cultural field  at  the  M2  site  (fig.  1). 
Samples  collected  from  wells  perforated 
at  10  and  6  m  below  land  surface  are  the 
most  evaporated.  They  also  are  the  most 
saline  and  have  the  highest  selenium  con- 
centrations (median,  3,500  yg/L) .  Sam- 
ples from  wells  at  the  3-m  depth  are  less 
evaporated  and  have  lower  selenium  con- 
centrations (median,  1,600  ug/L) .  This 
water  has  displaced  the  enriched  water 
since  the  mid-1970' s,  when  this  area  had 
a  water  table  within  1  m  of  land  surface. 
The  present-day  water  table  fluctuates 
between  3  and  1.5  m  of  land  surface. 

The  points  plotted  in  figure  10 
represent  an  evaporative  trend  that  is 
described  by  the  equation: 


6D  =  2.80x6190-40.7  (r=0.92), 


(8) 


The  slope  of  this  line  is  significantly 
(a=0.01)  lower  than  the  evaporative  trend 
slope  of  4.28  reported  by  Deverel  and 


Fujii  (1988)  for  shallow  ground  water 
collected  beneath  agricultural  fields  in 
the  western  valley,  north  of  the  M2  site. 
Allison  (1982)  presented  evidence  that 
indicates  that  this  low  slope  is  the  re- 
sult of  prolonged  evaporation  from  a  dry 
surface.  The  slope  for  six  unsaturated- 
zone  isotope  values  collected  from  this 
field  is  2.38.  This  and  evidence  pre- 
sented by  Allison  (1982)  indicates  that 
displacement  since  the  mid-1970 's  of 
unsaturated-zone  water  that  has  been  sub- 
ject to  evaporation  from  a  dry  surface 
soil  resulted  in  the  present-day  pattern 
of  isotopic  enrichment  in  the  shallow 
ground  water  shown  in  figure  10. 

Samples  collected  from  wells  at  a  depth 
of  12  m  are  similar  in  isotopic  composi- 
tion to  the  samples  collected  in  wells  at 
a  depth  of  3  m.  Because  the  isotopic 
composition  of  this  water  is  similar  to 
the  irrigation  water  applied  prior  to 
1970,  the  water  at  12  m  is  probably  water 
that  was  applied  at  land  surface  25  to  40 
years  ago.  Tritium  data  indicates  that 
vertical  rates  of  ground-water  movement 
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in  this  field  have  been  0.3  to  0.75  m/a 
during  the  last  20  years.  The  high  sele- 
nium concentrations  in  this  ground  water 
are  probably  the  result  of  leaching  of 
soluble  seleniviin  from  the  unsaturated 
zone  when  this  area  was  initially  irri- 
gated. The  local  well  water  has  less 
than  10  ijg/L  of  selenium.  Tritium  is  not 
present  in  samples  collected  below  6  m 
indicating  that  this  water  infiltrated 
prior  to  1970.  A  rise  in  the  water  table 
caused  the  water  now  present  at  the  6  and 
10  m  depths  to  become  evaporated,  in- 
creasing its  isotopic  enrichment  and 
selenium  concentration.  Installation  of 
drains  in  1980  has  since  lowered  the 
water  table  so  that  samples  collected  at 
the  3-m  depth  (near  the  water-table 
surface)  are  less  evaporated. 


Mineral  Controls 

Trends  in  major-ion  composition  corre- 
spond to  differences  in  ground-water 
salinity  and  selenium  concentrations 
(fig.  11)  .  Sodium  and  sulfate  are  the 
predominant  ions  in  the  most  saline  water 
(those  having  dissolved  solids  greater 
than  6,500  mg/L) .  The  chemistry  of 
ground  water  that  is  less  saline  is  more 
dominated  by  calcium  and  bicarbonate. 
Dissolution  and  precipitation  of  gypsum 
and  calcite  have  probably  resulted  in 
less  saline  calcium  sulfate  and  calcium 
bicarbonate  water. 

All  the  water  samples  were  super- 
saturated with  calcite  as  determined  by 
WATEQF  simulation  (Plummer  et  al.,  1976). 
The  increase  in  bicarbonate  percentages 
with  decreasing  salinity  is  probably  a 
manifestation  of  calcite  equilibrium  not 
apparent  in  water  having  higher  salini- 
ties.  High  partial  pressures  of  carbon 


dioxide  in  the  root  zone  also  may  lead  to 
higher  bicarbonate  percentages  also  not 
apparent  at  higher  salinities. 

The  primary  mineral-water  interaction 
in  the  alluvial  shallow  aquifer  is  the 
dissolution  or  precipitation  of  gypsum. 
As  salinity  and  isotopic  enrichment 
increase,  gypsum  saturation  indices 
increase  (fig.  12)  .  In  the  ephemeral- 
stream  alluvial  fans  and  the  margins  of 
the  four  major  fans,  ground  water  is 
generally  saturated  or  supersaturated 
with  gypsum.  In  nonsaline  areas  of  the 
major  fans,  ground  water  is  generally 
undersaturated  with  gypsum.  Gypsum 
precipitation  and  dissolution  are  the 
major  processes  controlling  sulfate 
concentrations.  Lower  salinity  values 
are  associated  with  lower  sulfate  activi- 
ties. As  salinity  increases  because  of 
evaporative  concentration,  gypsum  precip- 
itation limits  the  increase  in  sulfate 
concentration . 

Both  sulfate  and  selenium  increase  as 
salinity  increases  but  selenium  concen- 
trations increase  disproportionally  to 
sulfate  concentrations.  In  saline, 
selenium-enriched  areas,  selenium- to- 
sulfate  ratios  are  higher  compared  with 
the  low-selenium  areas.  In  the  upper  and 
central  areas  of  the  alluvial  fans 
deposited  by  intermittent  streams, 
selenium-to-sulfate  ratios  are  about 
0.00001.  In  more  saline  and  high  sele- 
nium areas,  such  as  the  lower-fan  areas 
of  the  intermittent-stream  fans  and 
middle-  and  lower-fan  areas  of  the 
ephemeral-stream  fans,  ratios  range  from 
0.0001  to  0.001.  Similarly,  Deverel  and 
Fujii  (1988)  found  that  the  selenium-to- 
sulfate  ratio  increased  exponentially 
with  increasing  isotopic  enrichment  of 
the  water  sample. 
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FIGURE  11.     Chemical  composition  of  water  samples. 
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FIGURE  12.   Comparison  of  gypsum  saturation  indices 
and  dissolved  solids. 


In  water  with  low-selenium  concen- 
trations (less  than  200  pg/L)  anc3  salin- 
ities, sulfate  and  selenium  probably 
behave  similarly.  In  the  low  range,  both 
sulfate  and  selenium  probably  are  simi- 
larly affected  or  unaffected  by  mineral 
precipitation,  dissolution,  or  other 
reactions  with  mineral  surfaces.  Leckie 
et  al.  (1979)  and  Kingston  et  al.  (1974) 
demonstrated  that  selenium  in  the  sele- 
nate  (+5)  oxidation  state  adsorbs  weakly 
to  mineral  surfaces  such  as  smectite 
clays  and  iron  oxyhydroxide  surfaces.  As 
ground-water  salinity  increases  because 
of  evaporation,  selenium  is  less  affected 
than  sulfate  by  mineral  precipitation  or 
surface  phenomena,  whereas  increases  in 
sulfate  activities  probably  are  limited 
by  gypsum  precipitation. 


Oxidation  Reduction 

The  oxidation-reduction  status  of  the 
46  samples  collected  during  spring  1986 
is  shown  in  figure  13.   The  pH  and  two 
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FIGURE  13.  Oxidation  reduction  status  of  selected  samples. 
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redox  indicators',  dissolved  oxygen  and 
the  selenate/selenite  couple,  indicate 
oxidized  and  alkaline  conditions.  Fer- 
rous iron  was  not  detected  at  concentra- 
tions greater  than  10  yg/L,  which  further 
indicates  oxidized  conditions.  Ferric- 
iron  concentrations  ranged  from  60  to  100 
yg/L.  Dissolved-organic-carbon  concen- 
trations were  high  (Thurman,  1985),  with 
a  median  of  4.9  mg/L,  but  organic-matter 
content  of  the  alluvial  materials  gener- 
ally is  low  (usually  less  than  0.25  per- 
cent). The  origin  and  identity  of  the 
dissolved-organic  matter  is  currently 
being  investigated. 

The  platinum  electrode  indicates  more 
reducing  conditions  and  is  in  disequilib- 
rium with  the  selenium  species  in  solu- 
tion. This  may  be  indicative  of  more 
reducing  conditions  due  to  some  other 
redox  couple.  The  oxygen-water  couple  is 
typically  grossly  out  of  equilibrium 
(Ball  et  al.,  1979).  The  presence  of 
some  selenite  in  solution  may  indicate 
that  nonassimilatory  reduction  of  sele- 
nate  coupled  with  the  oxidation  of 
organic  carbon  or  another  inorganic  con- 
stituent is  taking  place.  The  activity 
of  the  selenite  species  in  solution  is 
low  because  it  is  rapidly  adsorbed  onto 
mineral  surfaces.  This  hypothesis  is 
partly  supported  by  data  presented  by 
Fujii  et  al.  (1988).  They  reported  that 
about  50  percent  of  the  selenium  adsorbed 
on  soil  samples  collected  in  the  study 
was  in  the  selenite  form. 


CONCLUSIONS 


The  results  of  this  study  illustrate 
the  processes  that  have  resulted  in  the 
regional  distribution  of  salinity  and 
selenium  in  shallow  ground  water.  The 
main  conclusions  are  as  follows: 


(1)  The  highest  salinity  and  selenium  in 
shallow  ground  water  occur  in  alluvium 
deposited  by  ephemeral  streams  and  at 
the  margins  of  the  major  alluvial  fans, 
where  there  were  naturally  saline, 
fine-grained  soils. 

(2)  Low-to-moderate  salinity  and  selenium 
in  shallow  ground  water  are  associated 
with  upper  and  middle  areas  of  the 
major  alluvial  fans  deposited  by 
intermittent  streams. 

(3)  The  history  of  irrigation  in  the 
western  valley  has  affected  the  distri- 
bution of  soil  and  ground-water  salin- 
ity and  selenium.  Many  areas  with  the 
most  saline  soils  have  been  irrigated 
only  in  the  last  40  years.  These  are 
now  the  areas  of  highest  dissolved- 
solids  and  selenium  concentrations  in 
shallow  ground  water. 

(4)  In  the  area  irrigated  since  1940, 
present-day  ground-water  salinity  is 
significantly  correlated  with  natural 
soil  salinity. 

(5)  Oxygen-  and  hydrogen-isotope  data 
show  that  the  most  saline  and  highest 
selenium  water  present  at  low  altitudes 
resulted  from  evaporative  concentration 
when  the  water  table  was  near  the  land 
surface. 

(6)  Saline  and  selenium-enriched  water  is 
predominately  sodium  sulfate.  Water 
that  is  less  saline  is  calcium  sulfate 
and  calcium  bicarbonate  in  composition. 

(7)  Most  dissolved  selenium  in  ground 
water  is  in  the  selenate  form  and  the 
shallow  ground-water  samples  generally 
are  oxic. 
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